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ABSTRACT 


This  oociunent  has  been  approved 
Icr  public  release  and  sale;  its 
distribution  is  unlimited. 


With  intense  relativistic  electron  beams,  such  as  generated  by  induction 
linear  accelerators,  pulsed  transmission  line  accelerators  or  high  current  RF 

Cb~linacs,  collective  effects  can  play  an  important  role  in  the  FEL  process. 

CL_ 

^  High  gain  FELs  are  said  to  operate  in  the  Ran an  regime  when  collective  effects 


dominate  the  usual  ponderomotive  effects.  Collective  effects  may  also  play  at. 

; — I 

important  role  in  low  gain  FELs  when  the  bear-  current  is  sufficiently  high. 
O  In  this  paper  Lhe  linear  and  nonlinear  dynamics  of  the  FF.L  mechanisms, 

CS5  including  space  charge  effects,  is  analyzed.  The  various  collective  and 
noncollective  operating  regimes  are  discussed. 
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I.  Introduction 


Although  the  first  ideas  for  the  free  electron  laser  have  been  known 
since  the  original  article  by  Motz  [1]  in  1951  and  the  successful  experiments 
by  Phillips  [2]  in  1960,  it  is  the  recent  work  that  has  made  the  free  electron 
laser  a  serious  candidate  for  a  powerful  new  radiation  source.  In  fact  the 
free  electron  laser  has  become  the  conceptual  alternative  for  virtually  every 
radiation  source  from  the  microwave  to  the  ultraviolet  [3],  This  can  be 
understood  by  examining  Fig.  1,  where  the  power  vs.  wavelength  of  high  power 
coherent  radiation  sources  is  plotted.  These  sources  include  gas  and  solid 
state  lasers,  conventional  microwave  tubes  and  gyrotrons.  The  potential 
operating  range  not  only  covers  the  entire  wavelength  range,  but  maximum  power 
levels  are  comparable  with  the  other  conventional  high  power  coherent 
radiation  sources. 

Existing  accelerator  technology  together  with  the  characteristics  of  the 
FEL  Interaction  mechanism  divide  FELs  into  distinct  categories  [4-7],  '..'hese 
categories  are  distinguished  primarily  by  the  type  and  characteristics  of  the 
electron  beam  source. 

Free  electron  lasers  based  on  such  beams  as  RF  linacs,  microtrons  or 
storage  rings  can  be  expected  to  operate  in  what  Is  referred  to  as  the 


"Compton"  regime  [4-11].  Such  beams  are  generally  of  high  energy,  low  current 
and  high  quality  (low  emittance).  The  Compton  regime  is  one  in  which  the 
interaction  physics  is  primarily  governed  by  single-particle  effects; 
collective  or  space  charge  effects  can  usually  be  neglected.  Typically  In 
this  regime  the  radiation  gain  is  low,  thus,  practical  sources  operating  in 
this  regime  would  necessarily  function  as  oscillators  where  high  gain  is  not  a 
crucial  requirement.  In  the  absence  of  efficiency  enhancement  techniques,  the 
oper-  ng  efficiencies  are  generally  low,  (e.g.,  a  fraction  of  a  percent). 
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Since  the  beam  energy  and  quality  is  generally  high,  FELs  in  this  regime  can 
operate  in  the  optical  regime  or  beyond. 


Free  electron  lasers  based  on  intense  relativistic  electron  beams  (IREB) 
such  as,  pulse  line  accelerators  [12-15]  or  induction  linac  accelerators  [16- 
18],  operate  in  the  "Raman"  or  "collective"  regime  [4-7,19,20].  Here 
collective  effects  play  an  important  role  in  determining  such  characteristics 
as  the  radiation  growth  rate,  interaction  efficiency,  etc.  The  FEL  operating 
wavelength,  however,  remains  well-approximated  by  the  usual  expression 
appropriate  for  the  Compton  regime.  Numerous  FEL  experiments  have  been 
performed  with  pulse  line  generated  beams.  These  beams  are  produced  from 
plasma  induced  field  emission  diodes,  and  have  a  relatively  flat  voltage  and 
current  pulse  lasting  for  a  few  tens  of  nanoseconds.  Typically  they  are  in 
the  MeV  energy  range  and  carry  kiloamperes  of  current.  The  low  energy  and 
quality  of  these  beams  limit  their  operation  in  FELs  to  the  millimeter 
regime.  Since  the  beam  current  is  high,  the  radiation  gain  (or  spatial  growth 
rate)  can  be  large  enough  to  make  operation  as  an  FEL  amplifier  possible. 

There  Is  a  third  operating  regime  which  has  features  that  are  common  to 
both  the  Compton  and  Raman  regimes.  Here  the  wiggler  field  is  strong  enough 
so  that  the  ponderomotive  force  on  the  electrons  completely  dominates  the 
space  charge  forces  and  the  radiation  growth  rate  is  large.  We  will  refer  to 
this  regime  as  the  "High-Gain  Compton"  [4,5].  Free  electron  lasers  using 
storage  ring  generated  electron  beams  m3y  operate  in  this  regime  to  produce  x- 
rays  [21-23] . 

II  Linear  Theory  of  the  Free  Electron  Laser 

The  essential  feature  of  the  FEL  mechanism  is  that  the  electrons  undergo 
axial  bunching  in  the  combined  wiggler  and  radiation  fields.  In  this  section 


we  wilL  discuss  the  details  of  the  FEL  mechanism  in  the  linear  or  small  signal 
regime.  The  nonlinear  dynamics  will  be  discussed  in  section  III. 

The  physical  model  consists  of  a  relativistic  electron  beam  of  arbitary 
intensity  entering  and  propagating  through  a  static  helically  polarized 
wiggler  field.  Only  spatial  variations  along  the  z  axis  will  be  considered 
for  the  electron  beam,  wiggler  field,  radiation  fields  and  space  charge 
waves.  Figure  2  shows  the  basic  components  of  a  FEL  employing  a  linearly 
polarized  wiggler  field.  In  both  the  linear  and  nonlinear  analysis  contained 
in  this  paper  a  helically  polarized  wiggler  field  is  used.  Using  our  model  we 
obtain  expressions  for  the  radiation  field,  the  perturbed  charge  density  and 
ponderomotive  potential.  This  set  of  coupled  equations  fully  describe  the  FEL 
interaction.  A  dispersion  relation,  in  the  high  gain  regimes,  Is  derived 
which  shows  explicitly  the  coupling  of  the  electromagnetic  and  space  charge 
waves  by  Lhe  wiggler  field. 

For  cold  electron  beams,  we  find  two  regimes  where  the  radiation  grows 
exponentially.  In  the  collective,  or  Raman  regime,  the  contribution  to  the 
perturbed  beam  density  (and  hence  radiation  output)  due  to  the  self  consistent 
space  charge  potential  is  greater  than  that  due  to  the  ponderomotive 
potential.  In  the  high  gain  Compton  regime  (strong  wiggler  and/or  tenuous 
beam  limit)  the  contribution  to  the  perturbed  beam  density  from  Lhe 
ponderomotive  wave  dominates  that  due  to  the  self  consistent  space  charge 
potential.  This  regime  can  be  obtained  by  Increasing  the  wiggler  amplitude  or 
reducing  the  beam  density.  The  interaction  In  the  high  gain  Compton  regime  Is 
not  a  collective  one  as  it  is  in  the  Raman  regime. 

We  also  obtain  Lhe  small  signal  radiation  gain  with  space  charge  effects 
in  the  low  gain  Compton  regime.  This  result  Is  obtained  by  taking  Laplace 
Iranforms  of  Lhe  three  equations  for  Lhe  ponderomotive  potential,  density 


perturbation  and  radiation  field.  The  low  gain  Compton  mechanism  is  a  result 
of  constructive  interference  of  modes  rather  than  exponential  growth  of  an 
instability. 

The  helically  polarized,  static,  periodic,  magnetic  wiggler  is  taken  to 
have  the  form 


B^z)  =  Bw  (cos(kwz)ex  +  sin  (kwz)ey), 


(1) 


where  Bw  is  constant,  kw  -  2ir /Xw,  and  X  is  the  wiggler  wavelength.  This  form 

is  sufficient  for  particles  near  the  axis,  i.e.,  k  r,  «  1  where  r,  is  the 

w  b  b 

beam  radius.  A  more  realistic  representation  can  be  found  in  Refs.  8  and  9. 
The  vector  potential  associated  with  this  wiggler  field  is 

=  Aw(eikwz  e_  +  e  ~^wz  e+),  (2) 


where  Aw  »  Bw/kw,  and  e±  ■  (ex  ±  Electrons  streaming  axially  in  the 

presence  of  the  wiggler  field  (2)  and  radiation  field  produce  a  driving 
current  which  can  stimulate  (amplify)  the  imposed  radiation  field.  The 
radiation  field,  which  we  will  represent  by  its  vector  potential  A  ( z , t ) , 
evolves  according  to  the  wave  equation 


4 


(3) 


where  J_^  is  the  ponderomotive  induced  transverse  driving  current.  In  a  form 
similar  to  the  representation  of  the  wiggler  field  (2)  we  represent  the 


4 


stimulated  radiation  by 


(4) 


,(2,0  =  AR(eic(z>c)a+  +  e-^’^e..), 


where  A^  is  the  amplitude,  'Kz.t)  =  kz-u)t  is  the  phase,  k  is  the  complex 
wavenumber  and  as  is  the  real  frequency.  The  wiggler  induced  transverse 
driving  current  is  given  by 


J,(z,t)  =  -  |  e  |  (<5nv  +  n  vO  +  (r.onresonant  terms),  (5) 

where  <5n  is  the  perturbed  beam  density,  _v  is  the  transverse  velocity  induced 
by  the  wiggler  field,  nQ  is  the  ambient  beam  density  and  _v  is  the  transverse 
velocity  induced  by  the  radiation  field.  The  transverse  wiggle  and  radiation 
velocities  are  obtained  by  the  relativistic  force  equation  for  the  particles, 


dP/dt  -  -  S  e|  (E  +  (P  x  B)/YmQc),  (6) 

2  1/2  -1/2 

where  =  ym^y,  y  -  (1  +  P/P/n^c  )  =  (1  -  v*v/c)  is  the  total  relativ¬ 

istic  mass  factor.  The  fields  E  and  B  are 


E(z,t) 


3$  ~  _  i  »Ar 

"5z  ez~  c 


(7a) 


1  (  z ,  t ) 


(7b) 


where  <5  is  the  space  charge  potential  associated  with  the  perturbed  density, 
6n,  and  is  given  by 


From  (6),  conservation  of  canonical  transverse  momentum,  i.e., 
Ym o v  -  |  e |  A^/c  =  constant,  implies  that 


— w  i®  c  — « 

o 


and 


(8a) 


^R  yin  c  — R  * 


(8b) 


Substituting  (5)  together  with  (8)  into  (3)  yields 


2 

c 


*  4 


4ir  |  e  1 

2 

Y  m  c 
o  o 


A  . 
— w 


(9) 


From  Eq.  (9),  we  see  that  the  radiation  field,  A^,  is  driven  by  the  wiggler 
field,  A^,  and  the  beam  density  perturbation,  5n.  The  perturbed  beam  density 
must  now  be  determined  self  consistently.  From  charge  conservation,  the 
perturbed  beam  density  is  given  by 


3<5n 

3t 


W 


(10) 


where  ?J  is  the  perturbed  axial  beam  current  given  by 


6j  (z,t)  =  -  lei  (n  6v  +  5nv  ). 
z  ’  1  1  o  z  zo 


(11) 


In  Eq .  (II),  <$v^  and  v  are  the  perturbed  and  unperturbed  axial  electron 
velocities.  Combining  (10)  and  (11)  yields  the  following  expression  for  the 
perturbed  density, 
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Taking  the  axial  component  of  the  force  equation,  Eq .  (6),  and  using  the 

2 

relation  dy/dt  =  -  |e|(v*E)/m  c  we  find  that 

1  1  —  —  n 


dv  I!  ix  (v  x  B)»e 

_5.  =,  _  M  /-  11  +  ~  ~  2 

dt  ym  v  3z  c 

o 


v  (v  •  E) 

Z  —  -  , 


Linearizing  (13)  by  keeping  terms  to  first  order  in  the  radiation  field  yields 


--■J4 


In  Eq .  (14),  the  Yz  relativistic  reduction  in  the  space  charge  field  comes 

2 

from  combining  the  terms  3<5/3z  and  -(v  /c)  3$/3z  in  Eq .  (13)  where 

zo 

2  2  -1/2 

Y,  =  (1  -  v  /c  )  .  The  relativistic  factors  y  and  y  are  related  by 

z  z  z 

Y  =Y(l+(|e|A/mc)J  .  The  axial  force  terms  due  to  v  x  b/c  and 

2 

V  (v  ,•  E.)/c  can  be  written  in  terms  of  an  effective  potential  called  the 

zo  — l  —x  r 

ponderoraotive  potential,  <5^,  where 

,  ,  -Iel  V4  "le!  AwAR  i((k  +  k  )z-0»t) 

4  (z,t)  =  - =  - **  e  +  c.c,.  (15) 

P  Y  m  c  2y  m  c 

O  O  0  0 

We  see  from  Eq .  (14)  that  the  beam  velocity  and,  hence,  density  is  driven  by 
both  the  ponderoraotive  and  space  charge  potential  waves.  The  ponderoraotive 
wave  in  turn  is  proportional  to  the  wiggler  and  radiation  field  amplitudes. 
Taking  the  convective  time  derivative  of  both  sides  of  (12)  and  employing  (14) 
yields 


1_  (|_  +  |_  >*  (Zj0 


3z  3z 
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2  2 

Substituting  (7c)  into  (16)  and  using  3  £/3z  =  4tt  |  e  j  6 n  we  obtain 


d2Sn  ,  “p  .  ",e|no  3  ,3  Vzo  3  H  (17) 

7 T~  +  2^  5n  “  ~m  37  (3l  +  ~  3t)$pU’t)* 

dt  y  y  oo  c  y 

o  z 

Equation  (17)  shows  that  the  perturbed  charge  density  is  driven  by  the 
ponderomotive  potential  wave.  Equations  (9)  and  (17)  together  with  (15)  form 
a  set  of  coupled  relations  for  the  radiation  field  and  perturbed  charge 
density.  The  beam  density  perturbation  is  driven  by  the  ponderomotive  wave, 
which  is  proportional  to  the  radiation  field.  The  radiation  field  in  turn  is 
driven  by  both  the  ponderomotive  wave  and  wiggler  field.  The  coupling  in 
these  equations,  under  appropriate  conditions,  can  lead  to  radiation  growth. 


FEL  Dispersion  RelaLion 

The  phase  of  the  ponderomotive  wave  is  (k  +  kw)z~wt.  From  Eq.  (17)  we  see 
that  the  perturbed  density  should  have  a  similar  dependence  in  the  time 
asymptotic  limit,  hence  we  write 


x  ,  .  v  i((k  +  k  )z-u>t) 

9n(z,t)  =  6n(k,w)e  w  +  c.c, 


Using  (18)  together  with  (15,),  Eq.  (17)  becomes 


|(w-v  (k+k  ))2 - 2— r 

1  zo  w  L‘ 


w  A 

~f~2  "2  <k  +  V(k  +  kw 

8tty  m  c 
o  o 


Equations  (9)  and  (17)  can  be  used  in  the  high  gain  regime  to  obtain  a 


dispersion  relation  for  the  radiation  field.  In  this  section  transient  terns 


8 


associated  with  the  initial  value  solution  of  (9)  and  (19)  will  be  neglected. 
Substituting  (18)  together  with  the  representations  for  A^(z)  and  A^(z,t) 
given  in  (2)  and  (4),  into  (9)  gives 


2 


V 


—  4tt  e  | 


2 

■  fin  A  . 
w 


(20) 


Eliminating  fin  and  from  (19)  and  (20)  yields  the  following  dispersion 
relation 


l“2"  c2k2-^oJ  U“  -  vzo(k+kv»2  ■  “b/YoYzi 


2 


where  =  c6w  =  |  e  |  Ay/( y^c)  is  the  magnitude  of  the  electron  wiggle 

velocity.  In  obtaining  the  final  dispersion  relation  In  (21)  we  used  the 

2 

approximations  k  «  (1  +  0Z>YZ  kw  »  k  and  a>  =>  ck  to  simplify  the  terms  on  the 
right  hand  side  of  (19),  i.e., 

V  (1) 

(k  +  k  )  (k  +  k  -  -£2-  )  -  2  kk  . 
w  w  2  w 

c 

In  (21)  the  first  term  in  brackets  represents  the  uncoupled  electromagnetic 
node  while  the  second  bracketed  term  represents  the  two  uncoupled  beam  space 
charge  modes  having  an  effective  wavenumber  k  +  k  .  The  wiggler  field  provides 
for  the  coupling  between  the  electromagnetic  and  space  charge  modes.  Since  we 
will  be  primarily  concerned  with  the  forward  travelling  radiation  field  we  can 
approximate  the  electromagnetic  mode,  in  (21)  by 


9 


(22) 


2  2 
v  Y  Y 
zo  o  z 


2222  222  i  /"> 

n  -c  k  -uT  /y  “  2  u)(oj-  (c  k  +  u./y  .  '  "). 
bo  bo 


Using  (22)  the  FEL  dispersion  relation  becomes 


(  w2/c2  -  o^/yoc2)1/2)  [(k  +  kw  -  “/v20)2 


2  2 

VC  «2k“ 

^7  e»7" 

zo 


The  dispersion  relation  in  (23)  may  now  be  re-cast  into  the  form 


(k  -  k  ) (k  -  k  )(k  -  k  )  =  -  a  , 
em  -  + 


2  2 

2  \/c  7  kw 

where  a  =>  — 2 -  0  — —  is  the  coupling  coefficient, 

2To  W  B2 
zo 


k  -  (m2  -  w2/y  )1/2/c, 
em  bo 


(25a) 


is  the  electromagnetic  mode  wavenumber  and 


k+  3  vzokw  T  “b^VV^zo’ 


(25b) 


is  the  wavenumber  of  the  positive  and  negative  energy  beam  space  charge  mode 
respectively.  We  can  now  distinguish  two  operating  regimes  of  the  FEL 
mechanism,  often  referred  to  as  the  high  gain  r. encollective  (strong  wiggler) 


and  Raman  regimes. 
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High  Gain  Compton  Regime 

In  this  regime  the  forces  on  the  beam  electrons  due  to  the  ponderomot ive 
wave  dominates  that  due  to  the  collective  space  charge  effects.  This  can  also 
be  referred  to  as  the  strong  wiggler  field  limit  or  tenuous  beam  limit  and  is 
a  noncol lect ive  process.  In  this  limit  the  dispersion  relation  (23)  reduces 


(k  "  kem^k  ~  (a,/vzo  "  kw»2  =  "  “2‘ 


where  the  space  charge  term  which  comes  about  from  the  self  consistent  scalar 
potential,  i.e.,  «./( v  Y  /y )  has  been  neglected.  From  the  character  of  the 

U  U  o  u 

dispersion  relation  in  (26)  it  is  clear  that  this  regime  involves  the  coupling 

of  an  electromagnetic  mode  with  the  ponderomot ive  wave.  The  maximum  spatial 

growth  rate  occurs  when  k  is  equal  to  (&>  -  v  k  )/v  and  is  given  by 

em  zo  w  zo 


B2  w?  k  1/3 


r  /3  ,Bw  “b 

r  “  r  ^2 - T 

Y  c 
o 


In  obtaining  (27)  we  used  the  approximations,  co  »  ck  and  B  »  1.  The 

oz 

dispersion  relation  in  (26)  was  obtained  by  neglecting  the  space  charge  term, 
this  Implies  the  following  Inequality  between  the  wiggler  field  amplitude  and 
beam  density 


B  »  8  . 

w  cri 


,  b  .1/2 

t  C  1/2,  3 

cy  k  y 
o  w  z 


At  the  maximum  growth  rate  the  frequency  of  the  electromagnetic  wave  is 

2  2  _ 

w  ~  U  +  -zo)  Yz  vzo  kw  a  2Yz  ck^,  where  we  have  neglected  m^//yo  compared 


lO  '1). 


Raman  Regime 


The  Raman  regime  involves  a  coupling  between  the  forward  travelling 

electromagnetic  mode  and  the  negative  energy  beam  space  charge  mode.  In  this 

regime  the  beam  plasma  frequency  is  sufficiently  high  that  the  coupling 

between  the  electromagnetic  wave  and  the  Lwo  beam  waves,  i.e.,  negative  and 

positive  energy  modes,  can  be  considered  independently.  The  dispersion 

relation  describing  the  interaction  between  the  negative  energy  beam  wave  and 

electromagnetic  wave  is  obtained  from  (24).  Here,  the  effect  of  the  positive 

energy  beam  mode  on  the  coupling  is  weak  and  therefore  k  -  k+  can  be  replaced 

bv  (k  ~  k.)  =  2o), /(y  v  /y~ )  in  (24).  The  resulting  dispersion  relation  is 
J-'rbzzoo 


(k  -  k  )(k  -  k  )  -  - 
era  - 


cx^y  /y  v 

2  O  7.0 

2o>. 


(29) 


The  maximum  Raman  growth  rate  occurs  when  k  *  k  an(j  is 

sin  ■*  " 

uj  k  . 

r  =  3  (-h--j_'w) 1/2  .  (30) 

w  4/y  C 
o 

The  coupling  with  the  positive  energy  beam  wave  is  indeed  weak,  permitting  the 
replacement  of  k  -  k+  by  (k_  -  k+)  if 


e  «  (- 


2  m. 


w 


c/y  k  yJ 
o  w  z 


1/2 


(31) 


Note  that  (31)  is  the  reverse  inequality  as  in  the  High  Cain  Compton  Regime 
given  by  (28). 

Another  method  of  separating  the  Raman  regime  from  the  high  gain  Compton 
regime  is  that 
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w 

p 


»  CO 


p,crit 


,  1  /  4 


k 

w 


Y 


z 


in  the  Raman  regime,  while  co 

P 


co 

P.crit 


in  the  Compton  regime. 


Low  Gain  Compton  Regime 

In  this  section  we  will  simply  state  the  low  gain  expression  for  what  is 
usually  referred  to  as  the  Compton  regime.  In  this  limit  the  space  charge 
potential  usually  plays  a  small  role.  It  differs  from  the  high  gain  Compton 
regime  by  the  role  of  the  initial  conditions.  In  the  previous  section  we  took 
Fourier  transforms  to  obtain  the  FEL  dispersion  relation  in  (21).  The  Raman 
and  high  gain  Compton  regimes  exist  in  the  long  time  asymptotic  limit  and 
their  characteristics  are  independent  of  the  initial  conditions.  To  obtain 
the  gain  in  the  present  regime  we  take  Laplace  transforms  of  Eqs.  (9)  and 
(17).  The  low  gain  Compton  regime  gain  with  space  charge  effects  [24]  can  be 
shown  to  be  given  by 


G(z) 


8  2 
V 


k  z 
w 


3  le  (l  + 


„  2 
,sinB 

K  0  ;  ’ 


(32) 


1  /  2 

where  k,  =  (u>,/y  c)/(24y  )  ,  |G(z)l«  1,  0  *  Akz/2  and  Ak  =  k  +  k  -  u>/v  . 

b  b  zo  o  1  w  zo 

2  2 

In  the  absence  of  space  charge  effects,  i.e.,  k,  z  «  1,  the  small  signal 

b 

gain  in  (32)  reduces  to  the  usual  expression  [4-10,25].  Since  the  function 
2 

3(sin0/9)  / 3 0  has  a  minimum  value  of  0.54  when  9  =  1.3,  the  maximum  gain  is 
given  by 


G(z) 


max 


„2  2 

$  (D 

w  b 


w 


3 

z  > 


(33) 


where  z 


2.6/Ak  =  2 . 6/ ( k  +  k  -  tu/v  ) 


Estimates  of  the  laser  efficiency  for  monoenergetic  injected  electron 
beams  can  be  obtained  using  simple  trapping  arguments  [5,26],  In  the  FEL 
mechanism  electron  trapping  in  the  longitudinal  wave,  i.e.,  ponderomotive  and 
space  charge  wave,  is  the  saturation  mechanism  when  the  injected  electron  beam 
is  monoenergetic  in  the  longitudinal  direction.  The  axial  phase  velocity  of 
the  longitudinal  wave  is 


vph  =  «/(Re(k)  +  kw),  (34) 


where  Re(k)  is  determined  from  the  dispersion  relation.  In  the  linear 

development  of  the  laser  radiation  the  injected  axial  beam  velocity  is 

slightly  greater  than  the  phase  velocity,  v^  »  v  ^  +  Av  where 

v  »  Av  >  0  and  Av  depends  upon  the  particular  FEL  regime  under 
ph 

consideration.  Since  radiation  growth  occurs  when  Av  >  0,  the  phase  velocity 

•» 

of  the  longitudinal  wave  must  be  slightly  less  than  the  Initial  axial  electron 

velocity.  The  radiation  amplitude  increases  at  the  expense  of  the  electron's 

kinetic  energy  until  the  electrons  become  deeply  trapped  in  the  longitudinal 

wave.  At  this  point  the  radiation  field  reaches  its  maximum  amplitude  and  the 

average  axial  electron  velocity  is  approximately  given  by  vz  j  =  Vp^  -  Av. 

sat 

At  saturation  the  average  axial  electron  velocity  has  decreased  by 

approximately  2Av.  The  decrease  in  the  electron  beam  energy  can  be  directly 

equated  to  the  increase  in  radiation  energy.  For  highly  relativistic  electron 

2 

beams  the  decrease  in  the  average  electron  kinetic  energy  Is  Ae  =  2y  y  m  v  Av 

°  o  zo  o  zo 

ar.d  hence  the  radiation  efficiency  is 


n 


Ae 


(  r  -  1 )m  c 
o  o 


-  «  2y  Av/c. 
2  zo 


(35) 
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The  longitudinal  waves  "see"  the  beam  as  nonoenerget ic  if  the  beam's 

axial  velocity  spread  is  small  compared  to  Av.  Since  the  fractional  beam's 

2 

axial  energy  spread  is  Et^/EQ  =  Yzo  v^/c,  t'ie  :nonoenergetic  beam 
approximation  requires  that 


th  « 

1 —  <<  n* 
o 


(36) 


To  obtain  the  intrinsic  efficiency,  n,  at  the  maximum  growth  rate,  in  the  high 
gain  Compton  regime  we  solve  (26)  for  Re(k).  We  find 


w.  6  2/3  k 

Re(k)  =  to/v  -  k  +  ( - - — j 

02  w  77  ck 


w 


(37) 


Using  (37)  to  solve  for  Av  =  -  v  Eq.  (35)  yields  the  following 

expression  for  the  intrinsic  efficiency  in  the  high  gain  Compton  regime 


,  Vv  2/3 

n  =  (— - J 

4/7ockw 

! 

Following  the  same  procedure, 


(38) 

the  intrinsic  efficiency  in  the  Raman  regime  is 


n  = 


u>K//7 

D  O 

Y  ck 
z  w 


(39) 


and  in  the  low  gain  Compton  regime  is  given  by 

I 

n  =  X  / ( 2L)  (40) 

W 

I  where  L  is  the  interaction  length.  Table  I  lists  the  various  expressions  for 

the  spatial  growth  rates,  or  gain  and  the  corresponding  intrinsic  power 
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efficiencies  for  the  various  FEL  operating  regimes  discussed.  Operating  t 
FEL  at  shorter  wavelengths  by  increasing  the  beam  energy  or  decreasing  the 
wiggler  period  results  in  lower  efficiencies  and  more  stringent  requirements 
on  the  beam  energy  spread.  One  can  attempt  to  compensate  for  this  by 
increasing  the  beam  density,  however  this  is  usually  associated  with  an 
increase  in  beam  temperature. 


III.  Nonlinear  Theory  of  the  Free  Electron  Laser 

In  this  section  the  physical  model  is  identical  to  the  one  used  for  the 
linear  treatment  of  the  FEL,  except  here  the  wiggler  field  is  spatially 
tapered  [26-28].  Here  again  only  spatial  variations  along  the  z  axis  will  be 
considered  for  the  electron  beam,  wiggler,  radiation  and  space  charge  fields. 

The  variable  amplitude  and  period  wiggler  magnetic  field  can  be  expressed 
in  terms  of  the  vector  potential 


^(z)  =■  Aw(z){cos[j  kw(z')dz')ex  +  sinfj  kw(z')dz" )e  }  , 


(41) 


where  the  amplitude  A^z)  and  wavenumber  k^z)  are  known  and  slowly  varying 
functions  of  z.  The  potential  field  In  (41)  is  a  good  approximation  to  a 
right-handed  polarized  helical  magnetic  field  near  the  z-axis,  i.e.,  when 

<<  1.  where  r^  is  the  beam  radius.  The  wiggler  magnetic  field  associated 
with  (41)  is  given  by 


B  (z)  =  B  ( z ) { cos ( J  k  (z')dz'  +  a(z))e  +  sinfj  k  (z')dz'  +  a(z)|e  f,  (42) 

w  w  x  w  v 


where  3  (z)  =  -  ((k  (z)A  (z))2  +  (3A  (z)/3z)2 1 1  ^2,  and  a(z)  = 

w  W  w  W 


-l 


-  tan  J(3Aw(z)/3z)/(k  (z)A  (z))  ) ,  are  slowly  varying  functions  of  z.  The 


period  of  the  magnetic  field  is  X  (z)  =  2tt /( k  (z)  +  3a/3z)  ~  2tt /k  (z). 

w  w  w 

The  scattered  electromagnetic  and  electrostatic  fields  in  terms  of  the  vector 
potential  A^(z,t)  and  scalar  potential  $(z,t)  are  taken  to  be 

A^z.t)  =  Ar(z)  (cos  $(z,t)  ex  -  sin  d> ( z , t )  ey}  ,  (43a) 

$(z,t)  =  ®1(z)cos^p(z  ,t)  +  *2(z)sin<J>  (z,t),  (43b) 

z  z 

where  <t>(z,t)  =  J  k(z')dz'  -  at  +  0,  4>  (z,t)  =  J  (k(z')  +  k  (z')ldz'  -  at 

„  p  w 

o  o 

are  the  phases  and  the  amplitudes  of  the  potentials,  A^z),  ^(z)  and  $ 2(z) 
and  the  wavenumbers  k(z)  and  '^(z)  are  slowly  varying  functions  of  z.  The 
scattered  electromagnetic  field  represented  by  Eq .  (43a)  is  a  right-handed 
elliptically  polarized  field  traveling  in  the  positive  z  direction.  The 
frequency  a  of  the  field  and  the  phases  as  well  as  0  are  independent  of  z. 

Wave  Equations 

The  evolution  of  the  scattered  potentials  are  governed  by  the  wave 
equations 


2  2 

9 _ JL_  .2 — )A  (z  t)  = 

.  2  2  2JVZ’t) 

3z  c  3t 


4tt  .  ,  . 


and 


(44a) 


32*(z,L) 

”3*7 r~ 


4mJz(z,t)  , 


(44b) 


where  J(z,t)  is  the  driving  current  density.  Substituting  the  potentials  in 
Eqs.  (43)  into  (44)  yields 


(45a) 


2  2  2 

(u>  /c  -  k  (z)J  AR(z)  cos<j>( z , t ) 


-2k1/,2(z)  ~  (AR(z)kli/2(z)  j  sin$(z,t)  =  -  —  J  (z,t). 


•^■(^(z^incfi  (z,t)  -  $2(z)cos4>  (z,t))  =  ^-Jz(z,t), 


(45b) 


2  2 

where  terms  proportional  to  3  AR/3z  have  been  neglected.  The  coefficients  of 

the  sinusoidal  terms  on  the  left-hand  side  of  Eqs.  (45)  are  slowly  varying 

functions  of  z  and  independent  of  t.  The  arguments  of  the  sinusoidal  terms  on 

the  other  hand  are  rapidly  varying  function  of  t  for  z  fixed.  The  rapidly 

time-varying  terras,  in  for  example  Eq.  (45a),  can  be  removed  by  multiplying 

them  by  an<3  taking  the  temporal  average  over  one  wave  period, 

2ir/a) 

i.e.,  (u)/2ir)  J  dt.  Performing  this  operation  on  Eq .  (45a)  as  well  as 
o 

similar  operations  on  Eqs.  (45b)  yields 


’)')’)  ,  2it/u 

K/c  -  k  (z)  1  A^z)  =  -22  J  Jx(z,t)  cos*(z,t)dl, 

o 


(46a) 


2k^2(z)  -|j  (AR(z)k1/^2(z)  )  =  j  J^(z,t)  sin<t>(z,t)dt, 

o 


(46b) 


3$2(z) 


2u/u> 


(k(z)  +kw(z)j01(z) - ^ —  =  4]  Jz(z,t)  cos4>  (z,t)dt. 


(46c) 


3$.(z)  2ar  /  oj 

(k(z)  +  kw(z)  ]<>2(z)  +  — ^ —  =  4  J  Jz(z,t)  sin<f>  (z,t)dt. 

o  '  p 


3z 


(46d) 


Nonlinear  Driving  Currents 

It  is  now  necessary  to  obtaLn  expressions  for  the  transverse  and  axial 
components  of  the  current  densities  and  perform  the  time  integration  specified 
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in  Eqs.  (46).  In  general  the  non-thermal  electron  distribution  function, 
written  in  terms  of  the  electron  orbits,  is 


f(z,£,t)  =  nQVzo  j  6(z  -  z(to,t))S(px  -  px(tQ,t)) 


6(Py  -  Py^o’^)  6^PZ  "  Pz(to*t)ldt;o 


(47) 


where  nQ  is  the  uniform  particle  density  to  the  left  of  the  Interaction 

region,  i.e.,  z<0,  v  Is  the  constant  axial  electron  velocity  for 

zo 

z  <  0,  z(tQ,t)  is  the  axial  position  of  the  particle  at  time  t  which  crossed 
the  z  =  0  plane  at  time  tQ  and  £.(tQ,t)  is  the  momentum  vector  of  the  particle 
at  time  t  which  crossed  the  z  =  0  plane  at  time  tQ •  Thermal  effects  which  are 
characteristic  of  actual  electron  beams  can  be  easily  included  by  appropriately 
modifying  the  electron  distribution  function  in  (47).  The  integral  over  tQ 
in  Eq.  (47)  takes  into  account  the  continuous  flow  of  particles  Into  the 
interaction  region.  The  current  density  associated  with  this  electron 
distribution  Is 


J(z,t)  =  - - -  J  - - -  dt^  , 


m 


o  -°°  y(£(t  ,t)  )  j  3z(t  ,t)/3t| 


(48) 


where  y(j>)  =  (1  +  |.E^  1 /ro^c^)  *  ^  and 


dz" 


r(t  , z)  =  t  +  j  — r—~ 

O  O  '  V  (t  , 


o  z  o 


z')  ’ 


(49) 


is  the  time  it  takes  a  particle  to  reach  the  position  z  if  it  entered  the 
interaction  region,  z  =  0,  at  time  tQ  and  vz(tQ,z)  is  the  axial  velocity  of  a 
particle  at  position  z  which  was  at  z  =  0  at  time  t^. 


The  quantity  3z(t  ,t)/3t  is  the  axial  velocity  v  of  a  particle  at  time  t 
o  z 

which  was  at  z  =  0  at  Lime  tQ.  Clearly,  for  J_(z,t)  to  be  finite,  vz  should 
not  vanish  in  the  interaction  region.  We  assure  here  that  no  particle  is 
slowed  down  to  zero  velocity  in  the  laboratory  frame,  hence 
Y(£(to,t)  jmo|3z(tp,t^/at|  =*  PzCtQ,t)  and  the  driving  current  becomes 


J_(z,t)  =  - 


-  e  n  v 

o  zo 


06  2(t  ,t) 

J  FTTTO  5(t  "  T(to*2),dtc 


(50) 


Substituting  the  above  form  for  J(z,t)  into  the  right-hand  side  of  Eqs.  (46), 
we  obtain  the  self-consistent  amplitudes  and  phases  of  the  scattered 
potentials  in  terras  of  driving  currents. 

Since  the  system  of  particles  and  fields  are  in  the  temporal  steady- 
state,  particles  which  cross  the  z  =  0  plane  separated  in  time  by  2tt/w  will 
execute  identical  orbits  which  are  separated  in  time  by  2tt/w.  It  is, 
therefore,  possible  to  define  a  beam  segment,  "beamlet",  for  which  all 
possible  steady— state  orbits  of  the  actual  bean  particles  are  represented  by 
the  particles  In  the  beamlet,  but  are  displaced  in  time.  The  axial  length  of 
the  beamlet  is  clearly  Therefore,  substituting  (50)  Into  (46)  yields 


(w2/c2  -  k2(z)  )Ar( z ) 

v  2 tt/cu  p  J  t  ,x(t  ,z)  ) 

=  Meln0  —■  “  J  2 - ° -  cosd>fz,T(t  ,z))dt  ,  (51a) 

°  C  O  Pzfto,T(to,z))  °  0 


2k1/2(z)  —  (AR(z)k1/2(z) j 


V  2  IT  /  OJ  P  ft  ,  T  (  t  ,z)  1 

ZO  .  X  o  o 


=  -  4  |  e !  — —  w  J  — - -  sin.ifz,T(t  ,z))dt  , 

°  P.ft  *T(t  >z) )  °  0 


o  c 


(51b) 


z  o 
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f>(  z)  +  kv/z)  J^2)  ~ 


802(z)  2tt/u 

~a'z  =  "  4lelnoVzo  J  cosO  tz,T(Lo,z)ldto,  (51c) 

o  v 

30, (z)  2tr/aj 

(k(z)  +  kw(z)  ]02(z)  +  — ^ —  =  -  4|e|nQVzo  J  sln<i>  ( z, t ( tQ ,z) )dtQ  .  (  5 1  d) 

o  ” 

Notice  that  on  the  right-hand  side  of  the  above  equations  the  single  integrals 
over  tQ  are  from  0  to  2 n/tu.  As  we  will  illustrate,  these  integrals  can  be 
evaluated  numerically  by  following  the  orbits  of  a  relatively  small  number  of 
particles  which  enter  the  interaction  region  in  any  single  time  interval  of 
duration  2tt/o). 

Particle  Dynamics 

We  now  express  the  particle  orbits,  which  are  needed  for  the  evaluation 
of  Eqs .  (51)  in  terms  of  the  new  independent  variables  tQ  and  z.  The  forces 
exerted  on  the  electrons  arise  from  the  wiggle  and  potentials  given  in  Eqs. 
(41)  and  (43).  We  immediately  note  that  the  transverse  canonical  momenta  of 
the  particles  is  conserved.  Therefore,  if  both  the  wiggler  and  scattered 
fields  are  zero  as  z  *  -°°,  the  transverse  particle  momenta  are  given  by 

?x(2,t)  =  (A^z)  +  Ar( z , t )  J  •  ex,  (52a) 

and 

/ 

Py(z»t)  =  (A^z)  +  A^z.t))  •  ey.  (52b) 

Using  Eqs.  (52)  the  longitudinal  component  of  the  force  equation  can  be  put 


into  the  form 


2 


dp  (  z  ,  t  ) 
2 _ 

dt 


-  e 


2y(z ,t)mQC 


2  i|r:Vz)  + 


2 

-  2y(z ,t)  *(2,t)], 


(53) 


where  p2(z,t)  is  the  axial  momentum  and  the  relativistic  mass  factor  is 


2 

Y(z,t)  =  [1  +  (A^C z)  +  A^z.t))2  + 


ra  c 
o 


’2<Z't).l/2 

2  2  1 
m  c 
o 


(54) 


Equations  (52)  through  (54)  specify  the  particle  dynamics  in  terms  of  the 
wiggler  and  scattered  fields.  The  transverse  and  longitudinal  particle  motion 
is  decoupled.  Writing  Eqs.  (52)  and  (53)  in  terms  of  the  new  independent 
variables  z  and  tQ  we  find  that 


px(2,T)  „M  [A^cz)  ♦ 


(55a) 


(z,T)  =  (Atjv(z)  +  ARy(z,T)J, 


dpz(z,r) 

dz 


ll 

2c 


c  s  wy 

,2 


(55b) 


me  a 
O  ? 


2. — -t-ai-  !Vz)  +4(z’T)l  "  2y(z’t)-|7T"TzKz’t)1*  (55c) 

p  (z,t)  1  1 


We  have  expressed  the  particle  orbits  in  terms  of  the  entry  time  t  and  axial 

position  z.  Note  that  our  definition  of  the  momenta  implies  that, 

P  (t  »T)  =  Pv^z»T^»  P,/Lr,‘T)  =  P.,(Z»T)»  and  P„(t  »*)  =  p  (z,t).  To  obtain  the 

final  set  of  equations  for  the  amplitude  Ar(z)  and  wavenumber  k(z)  we  first 

combine  Eqs.  (51a)  and  (51b)  with  Eqs.  (51c)  and  (51d)  respectively.  Using 

•%*  ■>. 

the  expressions  for  p  and  p  given  by  Eqs.  (55a, b)  we  arrive  at  the  following 

x  y 

expressions 


,A.  W-h.nfc. 


a,: 
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(56a) 


-  ,  7  w  2ir/w 

(0)2/c  -  k  (z))Ar(2)  =J>_moVzo£  J  p“  ft0.T(t  0,Z)1 

4C  O 


{Aw(z)cos(4>  (z,T(to,z)l  +  A(z)}  dtQ. 


2k>'2(„ 

2c  o 


{A  (z)sin{<(>  (z,t(t  ,z)  J}dt  , 
w  p  o  o 


(56b) 


where  we  have  used  Eqs.  (41)  and  (43a)  for  A^(z)  and  AR(z,t)  and 

/(  .  .2  .  .1/2 

(*>b  =  (4n|e|  n^/nO  .  For  completeness  we  rewrite  Eqs.  (51c)  and  (51d)  for 
the  scalar  space  charge  potential 


3*2(z) 


(k(z)  +  kw(z)  J01(z) - ^ 


2  2 
-w,  v  me  2tt/u) 

T""?]7r  1  COSU  (z,T(to,z))dto, 
c  1  1  o 


(57a) 


3*,(z) 

(k(z)  +  kw(z))*2(z)  +  — j; — 


2  2 
-cu  v  me  2tt/oj 

”9 - 1  T  J  sin!  t>  f  z , t ( t  ,z)  !dt  . 

c2  11  |e|  jq  P  o’  o 


(57b) 


The  relevant  particle  dynamics  Is  contained  in  Eq .  (55b)  which  is 
rewritten  in  the  form 


dp  (t  , t)  i  [ 2  me..,. 

— JT —  =  4^“  lal  Uw(z)  +  Vz’x)1  "  2y(z’t)  17T~  ^“T^— J  <58> 

C  p  (t  t)  11 

Z  O 
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who  re 


P,(t  ,t)  1/2 


Y(Z,T)  =  [1  +  -IfLl^Cz)  +  A^z.-or  +  -~2~1 - j 


(59a) 


r(Ln‘z)  =  +  J  ~ 


z  y(z' ,r(t  ,z') J 


(59b) 


°  pz(to‘T(to*Z')) 


*ar<2) 


+  2A  (z)A(z)cosf  <p  (  z ,  t(  t  ,z))|. 
w  p  O' 


(59c) 


The  nonlinear  formulation  of  the  FEL  is  fully  described  by  Eqs.  (56),  through 
(58).  The  ponderomotive  potential  plays  a  central  role  in  axially  bunching 
the  electron.  From  Eq .  (58)  we  see  that  this  potential  is  given  by 


$pond(z’T)  =  ~2  Aw(z)AR(z)  cos(^  (z,x(to,z))!. 
Y  moc 


The  amplitude  and  phase  of  the  scattered  fields  as  well  as  the  axial  beam 
momentum  all  vary  with  a  characteristic  axial  length  which  is  much  longer  than 
the  viggler  wavelength  This  fact  allows  for  inexpensive  numerical 

simulations  to  be  performed  in  the  laboratory  frame  with  extremely  high  energy 
electron  beams. 


Coupled  Pendulum  and  Wave  Equations 

At  this  point  our  equations  can  be  simplified  and  'written  in  a  more 

*s» 

conventional  form.  To  this  end  we  define  the  electron's  phase,  i>,  with 
respect  to  the  ponderomotive  potential. 
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=  j 

o 


(k(z')  +  kw(Z')  -  U/V2(^0,z)’d2^ 


(61) 


where  =  ~  a)0t0  t*ie  e^ectron's  phase  upon  entering  the  interaction  region 

at  z  =  0  and  t  is  also  a  function  of  the  initial  value  3^/3 z  j  .  The  second 
derivative  of  the  electron's  phase  is  z=® 


d'X'l*  »z) 

7? 


f  (k 

dz  w 


+  k)  + 


0) 


v^(*  ,z) 

z  o 


dv  (\p  ,z) 
z  o 

dz 


(62) 


Substituting  (58)  into  (62)  and  assuming  y  »  1,  I A^ |  »  | A  |  and 
v^  a  u/(k  +  k  )  «  c  we  find  the  following  pendulum  equation  [29] 


djj)  d_ 


dz 


2  dd  +  k> 


I  1 2  k  ,  „ 

1§J _ w  1_  A  2 

2  4  2  dz  w 

moc  YX 


2  k‘ 


ym  c 


+  2  4  ~T  tAwAR  3in'*'  “  “Tif"  -  $2C0S^1  • 


(63) 


V  Y1 


f  2  1  /  2 

where  Y±  =  (l  +  (|e[Aw/mQc)  1  is  the  relativistic  mass  factor  associated 
with  the  transverse  wiggle  motion.  In  the  absence  of  space  charge  effects, 
(63)  reduces  to  the  usual  pendulum  equation  with  a  tapered  wiggler  [23]  and 
without  tapering  [10].  Using  the  same  assumptions  as  used  to  obtain  the 
pendulum  equation  we  find  that 


j  7  0  “>.  cos’4»(’)  ,z) 

/c  '  k‘(2))h<z)  ■  ~2  Vz>  <~!u7j  j  : 

c  o 


(64a) 


1/2,  x  d 
(z) 


dz  ‘R 


1  /  ° 

(7.W  U)  ]  = 


s  i  r. .  ■ 


2c 


f  AW(Z)  <“ 


z’0) 


(64b) 
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(64c) 


C 


r„ 


2  6.  w 

w,  m  c. 

Vz>  =  — 2  Tel-  <cos'^^0*z^>> 

a)  '  ' 

2  2 

2M,  m  c 

=  ~ — ifr  ^sin^'*'0»z^»  (64d) 

0)  I  I  c 

2ir 

where  <•••>  =  ]  («•«  )<j  /2ir  represents  the  ensemble  average  over  the 

o 

initial  phases  of  the  electrons.  Equations  (63)  and  (64)  represent  the  full 
set  of  nonlinear  equations  describing  the  one  dimensional  FEL  process  with 
space  charge  effects  and  tapered  wiggler. 

Numerical  results  in  the  Raman  regime 

In  Lhis  section  we  present  numerical  results  [27]  for  the  coupled 
nonlinear  FEL  equations  in  (63)  and  (64).  We  assume  that  a  monoenergetic 
electron  beam  enters  the  interaction  region  at  z  =  0  with  a  uniform  density. 
The  magnetic  wiggler  field  given  in  (41)  is  assumed  to  be  built  up 
adiabaticaliy  from  z  <  0  to  its  Initial  value  at  z  =  0.  In  all  of  our 
numerical  simulations  a  small  amplitude  radiation  field  is  introduced  as  a 
perturbation  at  z  =  0  and  allowed  to  grow  spatially  and  self-conslstently 
according  to  the  FEL  equations. 

We  first  consider  the  case  where  the  magnetic  wiggler  parameters  are 
fixed,  l.e.,  constant  amplitude  and  period.  Later  we  consider  the  case  where 
the  wiggler  period  is  adiabaticaliy  decreased,  resulting  in  substantially 
higher  radiation  efficiency. 

In  this  example  submillimeter  radiation  at  X  =  338  pm  is  generated  using 

a  2.6  MeV  electron  beam.  Table  II  lists  the  salienL  parameters  for  the 

wiggler  field  electron  beam  and  output  radiation.  The  magnetic  wiggler 

amplitude  2.5  kG  and  the  period  is  fixed  at  2.0  cm.  The  2.6  MeV  (y  =6),  5 

o 

k\  electron  beam  has  a  transverse  equilibrium  velocity  of  vw  =  0.078c. 
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Figure  5  shows  the  amplitude  of  the  vector  potential  of  the  excited 


radiation,  AR(z),  and  the  spatial  growth  rate,  r  =  JtnAR(z) 1 /8z,  as  a  function 

2  12  -1 

of  z.  These  plots  are  for  the  frequency  u>  =  2y  ck  =  5  x  10  sec 

oz  w 

(y  =  338  pm).  The  long  spatial  region  where  the  growth  rate  is  fairly 
constant  is  the  linear  region  of  the  interaction.  The  value  of  the  radiation 
frequency  in  this  figure  has  been  chosen  to  maximize  the  linear  growth  rate. 
Figure  6  shows  the  amplitude  of  the  beam  space  charge  wave  and  the 
ponderomotive  wave  as  a  function  of  the  interaction  distance.  Both  wave 

amplitudes  are  of  the  same  order  of  magnitude  signifying  that  the  process  is 

for  the  most  part  in  the  Raman  regime.  Figure  7  shows  the  evolution  of  the 
wavenumber  associated  with  the  radiation  field.  Since  the  radiation 

wavenumber  is  greater  than  m/c,  the  effective  index  of  refraction  in  the  beam 

region  is  greater  than  unity.  This  implies  that  the  radiation  field  will  tend 
to  focus  inward  towards  the  electron  beam.  Figure  8  shows  a  comparison 
between  the  linear  growth  rata  obtained  from  the  dispersion  relation  in  (24) 
(solid  curve)  with  the  growth  rate  obtained  from  the  linear  regime  of  the 
fully  nonlinear  simulations  (crosses  (x)).  Also  in  this  figure  the 
theoretical  efficiency  based  on  Eq.  (35)  (dashed  curve)  is  compared  with  the 
nonlinear  results  (circles  (o)). 


Efficiency  Enhancement 

The  phase  velocity  of  the  total  longitudinal  wave  potential,  i.e.. 


ponderomotive  plus  space  charge  is  approximately  v  ^  =  w/(k  +  kw) . 


The  longitudinal  wave  potential  is  responsible  for  axially  bunching  and 
eventually  trapping  the  electrons.  If  the  wiggler  period  is  held  fixed,  the 
radiation  field  reaches  its  maximum  value  when  the  electrons  are  trapped  in 
the  longitudinal  potential  wells.  Just  prior  to  the  saturation  of  the 


radiation  field,  the  electrons  are  somewhat  spatialLy  bunched  and  trapped  near 
the  bottom  of  the  wave  potential.  The  trapped  electrons  at  this  point  can  be 
considered,  for  our  purpose,  to  form  a  macro-particle.  By  appropriately 
reducing  the  phase  velocity  as  a  function  of  axial  distance  down  the 
interaction  region,  the  kinetic  energy  of  this  macro-particle  can  be  further 
reduced  and  converted  into  wave  energy.  The  phase  velocity  must  be  reduced  in 
such  a  way  so  that  the  inertial  potential  of  the  trapped  macro-particle  is 
always  less  than  the  potential  of  the  growing  longitudinal  wave.  The  phase 
velocity  can  be  reduced  by  decreasing  the  period  of  the  magnetic  pump  as  a 
function  of  z.  In  order  for  the  macro-particle  to  remain  trapped,  the  spatial 
rate  of  change  of  the  wiggler  period  must  be  sufficiently  slow.  In  principle 
virtually  all  the  kinetic  energy  of  the  macro-particle  can  be  extracted  and 
converted  to  wave  energy.  However,  not  all  the  beam  particles  comprise  the 
macro-particle;  some  are  untrapped.  Converting  particle  kinetic  energy  into 
radiation  by  varying  the  wave  velocity  is  somewhat  analogous  to  the  reverse 
process  of  particle  acceleration  in  say  an  RF  linac.  Figure  9  illustrates  the 
effect  on  the  saturated  radiation  amplitude  using  a  tapered  wiggler  field.  In 
this  figure  the  wiggler  period  is  gradually  decreased  spatially  at  the  point 
where  the  electrons  are  deeply  trapped.  Once  the  electrons  are  deeply  trapped 
they  are  decelereated  by  decreasing  the  wiggler  period  [27,28].  Figure  9 
demonstrates  that  dramatic  enhancements  in  the  FEL  efficiency  can  be  achieved 
in  this  way.  (The  parameters  of  Fig.  9  are  different  from  Table  II,  and  they 
are  given  in  Ref.  [27].)  There  are  of  course  other  schemes  for  enhancing 
efficiency  in  the  FEL.  These  include,  for  example,  spatially  decreasing  the 
wiggler  amplitude  [26,28]  and/or  applying  an  axial  accelerating  force  such  as 
a  D.C.  electric  field  [30,31]  to  the  trapped  electrons. 


As  a  concluding  comment,  concerning  the  rapid  growth  of  the  free  electron 
laser  concept,  the  theory  has  been  extended  to  include  axial  magnetic  fields 
[32,33]  thermal  effects  [34],  three-dimensional  effects  [35,36]  and  the 
nonlinear  stability  of  the  radiation  field  [37],  Also  a  number  of  FEL 
experiments  are  now  underway  or  in  the  planning  phase  and  numerous  workshops 
and  special  journal  issues  on  the  subject  have  appeared  [7,38-42]. 
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Table  I 


The  quantities  used  in  Table  I  have  the  following  definitions:  v  =  1/17  is 

Budker's  parameter,  I  is  the  beam  current  in  kiloanperes,  L  is  the  wiggler  length, 

is  the  beam  radius,  $w  =  v^/c,  YQZ  =  Yq/(1  +  ,f  (9)  =  3  (sin  9/0)2/30, 

0  =  (1  -  v  /v  )toj/2,  t  =  L/v  is  the  electron's  transit  time  and  F  is  the 
oz  ph  oz 

filling  factor,  i.e.,  beam  area/radiation  area,  (c.g.s.  units  are  used  unless 
otherwise  stated) 

FEL  Operating  Gain  or  Intrinsic  Power 

Regimes  Growth  Rate  Efficiency 


High  Galn- 
Compton 
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Table  II 


Raman  FEL  in  the  Submillimeter  Regime 
(Constant  Magnetic  Wiggler  Parameters) 


Magnetic 

Wiggler  Parameters 

Wiggler  Wavelength 

X 

w 

2.0  cm 

Wiggler  Amplitude 

Bw 

2.5  kG 

Electron  Beam  Parameters 

Beam  Energy 

Eo 

2.6  MeV(y  =  6) 
o 

Beam  Current 

xb 

5  kA 

Axial  Gamma 

Yzo 

5.4 

Beam  Radius 

rb 

0.3  cm 

Wiggle  Velocity 

8w 

0 .078 

Self  Potential  Energy  Spread 

aE/E0 

4.9  % 

Output  Radiation  Parameters 

Radiation  Wavelength 

X 

338  pm 

Linear  e-folding  length* 

l  =  i/r 

5.3  cm 

e 

Efficiency* 

n 

9.2  % 

Radiation  Power* 

P 

1.2  GIV 

*  At  the  maximum  growth  rate. 
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Figure  Captions 
Power  vs.  Wavelength 

The  power  vs.  wavelength  of  some  representative  high  power 
coherent  radiation  sources  are  shown.  The  output  from  recently 
reported  free  electron  laser  sources  are  shown  and  the  type  of 
electron  beam  source  used  indicated.  The  visible  FEL's  utilizing 
storage  ring  electron  beams,  are  at  a  power  level  below  the  scale 
of  the  graph. 

Shows  a  typical  FEL  configuration  employing  a  linearly  polarized 
wiggler  field 

Dispersion  diagram  of  the  coupled  beam  space  charge  modes  and 
electromagnetic  mode.  The  Raman  FEL  Instability  occurs  near  the 
intersection  between  the  negative  energy  beam  mode  and 
electromagnetic  mode. 

Spatial  growth  rate  and  intrinsic  efficiency  in  the  Raman  and  high 
gain  Compton  regimes  as  a  function  of  the  wiggle  veloctty/c. 

Shows  the  evolution  of  the  radiation  field  amplitude  and  the 
growth  rate  as  a  function  of  the  interaction  distance. 

Shows  the  spatiaL  evolution  of  the  beam  space  charge  wave  and 
ponderomot ive  wave  amplitude.  This  illustration  is  In  tha  Raman 


regime . 


Shows  the  spatial  evolution  of  the  radiation  wavenumber,  the  range 
where  k  is  constant  is  the  linear  regime.  Since  k  >  oj/c  the 
radiation  field  will  focus  inward  towards  the  electron  beam. 

Shows  the  comparison  between  the  linear  growth  rate  obtained  from 
the  dispersion  relation  (solid  curve)  and  that  obtained  from  the 
simulations  in  the  linear  regime  (crosses  (x)).  Also  shown  is  the 
theoretical  efficiency  based  on  Eq .  (35)  (dashed  curve)  compared 
with  the  nonlinear  simulations  (circles  (o)). 

Shows  the  enhancement  of  the  radiation  field  and  efficiency  by 
tapering  the  wlggler  period. 
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